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Abstract. We study the evolution of a A-body weakly interact- 
ing system of Bosons. Our work forms an extension of our pre- 
vious paper I [13] . in which we derived a second-order correction 
to a mean-field evolution law for coherent states in the presence 
of small interaction potential. Here, we remove the assumption of 
smallness of the interaction potential and prove global existence 
of solutions to the equation for the second-order correction. This 
implies an improved Fock-space estimate for our approximation of 
the A-body state. 



1. Introduction 

Experimental advances in the Bose-Einstein condensation (BEC) of 
dilute atomic gases [TJE] have stimulated interesting questions on the 
quantum theory of many-body systems. For broad reviews, see, e.g., 
[19l[2T]. In BEC, integer-spin atoms (Bosons) occupy macroscopically 
a quantum state (condensate). For a large number N of interacting 
atoms, the evolution of this system has been described fairly well by 
a single-particle nonlinear Schrodinger equation [151 HQ 1201 125] ■ The 
emergence of this mean-field description from the iV-body Hamiltonian 
evolution has been the subject of extensive studies; see, e.g., [5HT0 ] IT8 | 
[23]. 

In [13], henceforth referred to as paper I, we derived a new nonlin- 
ear Schrodinger equation that describes a second-order correction to a 
mean-field approximation for the A^-body Hamiltonian evolution. Our 
work was inspired by: (i) Fock-space estimates provided by Rodnianski 
and Schlein [23J, with regard to the rate of convergence for Hartree 
dynamics; and (ii) a second-order correction formulated by Wu [251126] . 
who introduced a kernel for the scattering of atoms in pairs from the 
condensate to other states. In paper I, we derived a new Fock-space 
estimate; and showed that for small interaction potential the equation 
for our second-order correction can be solved locally in time. 
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The present paper is a continuation of paper I. The main improve- 
ment presented here is the removal of our assumption on the smallness 
of the interaction potential. Notably, we prove global existence of so- 
lutions to the equation for the second-order correction. Our approach 
enables us to derive an improved with respect to time Fock-space esti- 
mate for our approximation of the iV-body quantum state. 

In the remainder of this introduction, we review elements of the 
Fock space, summarize the major results of paper I, and state the 
contributions of the present paper. For a more extensive discussion of 
the background, the reader may consult, e.g., the introduction in our 
paper I. 

Fock space and mean field. The problem at hand concerns the 
time evolution of N weakly interacting Bosons described by 

id t tp = H N tp , 

where ip is the iV-body wave function, the Hamiltonian operator 
H N : = J a* x A x a x dx - J v(x - y)a* x a*a x a y dxdy 

and v is the two-body interaction potential. A few comments on these 
expressions are in order. Here, we use the (convenient for our purposes) 
formalism of second quantization, where a*, a are annihilation and 
creation operators in a Fock space T [2], to be defined below; ip is a 
vector in J 7 ; and V is the particle interaction. Note that, in comparison 
to paper I, we changed the sign of the interaction term V, i.e., we 
replaced v with — v so that having v > corresponds to repulsive 
interaction, which leads to defocusing behavior. 

At this point, it is advisable to review the basics of the Fock space 
J 7 over L 2 (M. 3 ). For Bosons, the elements of T are vectors of the form 
V> = (ipo,ipi(xi), ip2(xi,x 2 ), ■■■ ), where ip <E C and ip n G L 2 s (R 3n ) are 
symmetric in x\, . . . , x n . The Hilbert space structure of T is given by 
the inner product ((/>, t/>) = J2 n I 'Pnipndx. 

For any / G L 2 (M 3 ), the (unbounded, closed, densely defined) cre- 
ation operator a*(f) : J 7 ' —> J 7 and annihilation operator a(f) : J 7 ' —> J 7 
are defined by 

(a*(/)^ n -i)(a;i, ■ ■ ■ , x n ) = — ^ f{xj)4>n-i{xi, Xj-i, x j+1 , ...,x n ) 

v " i 

7 = 1 



{a(f)iJ n+1 )(xi,x 2 , ...,x n ) = Vn + 1 / 4)(r l+ i)(x 1 x ll . . .,x n )f(x) dx . 
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The operator valued distributions a* and a x are defined by 

°*{f) = J f{ X ) a l dx , 

a(f) = J f(x) a x dx . 

It follows that the operators a, a* satisfy the commutation relations 

[a x , a* y ] = 5{x-y) , [a x , a y ) = [a*, a*] = . 

We are interested in the evolution of coherent states, i.e., vectors of 
the form e~^ A ^fl where Q = (1,0,...) G J 7 is the vacuum state, 
4>(t,x) is the one-particle wave function (to be determined later), and 



A(<f>) := j (0(x)a* — (f)(x)a x ) dx . (1) 
It is important to notice that 

£ -VNA(4>)Q =(...Cnf[ (j)(xj) ...) . 

Thus, the n th slot in the coherent state Fock vector consists of the tensor 

n 1 /2 

product of n functions 4>(x); the relevant constant is c n = (^f) 
Furthermore, the number operator, M := j a* x a x dx, satisfies 

(fte^W | M | e~^ A ^n) = NUW 2 . 

Thus, if we normalize the wave function by setting ||0|| = 1, the average 
number of particles remains constant, N. 

It can be claimed that a reasonable approximation for the many-body 
time evolution is expressed by the Fock vector 

./, — p VNAW))n 

Yappr ° ; 

where <j>(t,x) satisfies the Hartree equation ([3]). This ip ap p r encapsu- 
lates the mean field approximation for N weakly interacting Bosons. 
The precise meaning of this approximation as well as its rigorous jus- 
tification were studied within the PDE setting by Erdos, Schlein, Yau 
[pfflU] via Bogoliubov-Born-Green-Kirkwood-Yvon hierarchies for re- 
duced density matrices (see also Klainerman and Machedon, [18], for 
a simplification of the uniqueness part of the argument). In the Fock 
space setting, the mean field approximation was studied by Ginibre 
and Velo [12] and, most recently, by Rodnianski and Schlein [23]; see 
also Hepp, [T7j . 
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Main results of paper I. Starting with a coherent state as initial 
data, in [13J we proposed a correction of the form 

■- p -y/NA(4>) -B{k)Q 
Yappr ° ° 5 

where 

B(k) :— J (k(t, x,y)a* x a* — k(t, x, x)a x a y ) dxdy (2) 

and the kernel k(t,x,y) satisfies an appropriate nonlinear evolution 
equation. This k loosely corresponds to the "pair excitation function" 
introduced by Wu [25||26] but our set-up and derived equation for k are 
different. By assuming that the interaction potential v (x — y) is small, 
we proved that for a finite time interval our approximation stays close 
to the exact evolution in the Fock space norm. To be more precise, we 
proved the following general theorem. 

Theorem 1.1. Suppose that v is an even potential. Let <p be a smooth 
solution of the Hartree equation 

^ + A0-(^|0| 2 )0 = O (3) 

with initial conditions 0o- Assume all conditions (l)-(3) listed below: 

(1) The kernel k(t,x,y) G L 2 (dxdy) for all t, is symmetric, and 
solves the equation 

(iu t + ug T + gu - (1 + p)m) = (ip t + [g,p] + um)(l + p)~ 1 u , (4) 

where all products in (j4j) are interpreted as spatial compositions 
of kernels (or operator products), "1" is the identity operator, 
and 

u(t, x, y) := sh(k) := k + T^kkk + ■ • • j (5a) 

5{x — y) + p(t, x, y) := ch(/c) := 5{x — y) + ^kk + . . . , (5b) 

g(t, x, y) := -AJ(x - y) + v(x - y)<p(t, x)lj>(t, y) (5c) 

+ (v*\<f>\ 2 )(t,x)6(x-y) , (5d) 

m(t, x, y) := -v(x - y)~cj)(t, x)0(t, y) . (5e) 

(2) The functions defined by 

f{t) := \\e B \A,V}e- B ^iy , 

g(t) := He^e-^ll^ , 

are locally integrable; recall that V is the interaction operator, 
and A, B are operators defined by ([!]), (|2]). 
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(3) The trace f d(t, x, x) dx is locally integrable in time, where the 
kernel d(t, x, y) is 



d(t, x, y) = (ish(k) t + sh(k)g T + gsh(k)) sh(k) 
- (ich(k) t + [g, ch(k)]) ch(k) 




n 



Moreover, we showed that the hypotheses of this theorem are satis- 
fied locally in time if v is small. 

Theorem 1.2. Let to be sufficiently small and v(x) = x( x ) u?| f or X £ 
C^°(1R 3 ) . Assume that <f> is a smooth solution to the Hartree equation 
©, ||0||L2(<fa) = 1- Then, there exists k G L°°([0, l])L 2 (dxdy) that 
solves (jlj) for < t < 1 with initial condition k(0,x,y) = 0. In 
addition, we have the estimates 



Main results of this paper. In the present paper, we remove 
the smallness assumption on the interaction potential, prove that the 
evolution equation of k(t, x, y) has a global in time solution and obtain 
a stronger a priori estimate for the difference of the approximate and 
exact solution for the N-body Fock-space vector. In particular, we prove 
the following theorem. 

Theorem 1.3. Let the notation be as in Theorem \l.l\ Consider v(x) = 
> ; where x £ Co° and x(r) is a decreasing cut-off function. 

Assume 4>o has sufficiently many derivatives in L 2 and ||o;0o||l 2 — C. 
Further, suppose k(0, ■, ■) £ L 2 (M 6 ) is prescribed. Then, the hypotheses 
of Theorem \l.l\ are satisfied globally in time and 




and 





oll-F 



(7) 
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Remark 1.4. It follows from our calculations that if we omit the as- 
sumption v > 0, the hypotheses of Theorem 11.11 are still satisfied glob- 
ally in time, but we no longer have estimate ([7]). 

The remainder of this paper is largely devoted to the proof of Theo- 
rem [T2] and is organized as follows. In section [3] we derive the a priori 
estimate 



\u(T)\\ L 2 < (^j \\m\\ L 2dt+ |K0)|| L 2^ exp 




In section[2]we prove that J °° ||m||^2(it < C if v > 0. In section (J3J) we 
show that (jlj) is locally well posed for L 2 , possibly large, initial condi- 
tions for u. This proof is much harder than the corresponding one in 
paper I; the latter worked for zero (or small) L 2 initial conditions. The 
idea here is to transform the quasilinear equation (j4]) into an equivalent 
semilinear one. Section [5] is devoted to estimating the error terms / 
and g entering ([6]). In section [6] we construct the requisite operator e B 
in the case where ||/c||l 2 is large and e B is no longer defined as a conver- 
gent Taylor series; and elaborate on the connection of this construction 
with the Segal-Shale- Weil, or metaplectic, representation. Finally, the 
Appendix focuses on an improved computation of some error terms 
previously computed in section 8 of paper I. This leads to a simpler 
proof of our stronger estimate ([7]). Our notation is not uniform across 
sections, but is self-explanatory and convenient. When the variables 
are called x\ and X2, 4>i abbreviates 4>(xi), t>i_2 = v(xi — £2), etc. 



2. PSEUDOCONFORMAL TRANSFORMATION FOR HARTREE 

EQUATION 

The goal of this section is to find an estimate for the decay rate in 
time of \\4>(t, -)||l 4 (ir 3 ); where is a solution of the Hartree equation, 

^ + A0-(u*|0| 2 )0 = O, ( 8 ) 

with initial condition O such that H^oll// 1 + H^olU 2 be finite. For 
this purpose, we make use of the technology of dispersive estimates 
from [13] 

We start with some preliminaries. Let 



W = v* |0| 2 



(9) 
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The quantities relevant for the conservation laws (to be stated below) 
are defined by 

P-= (i/2)H 2 ; 

Pj := (l/2i) (0V,0 - 0V,0) ; po = (1/2*) (0<9 t - 0^0) ; 

Vjk ■= V j 0V fe + V fc 0Vj0 ; a j = Vj4>d t (f) + dtcjNjcj) . 

Let us call a := tr(ajk), the trace of the tensor cr^, and define two 
more quantities, namely, 

A : = —p + (l/2)<r + W> = Ap - Wp ; 
e : = (l/2)a + Wp . 

With regard to A, see (TTTT) . 

The quantity e is the energy density, while A is the Lagrangian den- 
sity. Indeed, one can see that the evolution equation can be derived as 
a variation of the integral 

£(0,0) : = J dx{\} . 

The associated conservation laws can be stated in the forms 

dtp - V,y = , (10a) 

dtPj - V fc {a/ - 5/A} + lj = , (10b) 

d t e - Vja j + l = . (10c) 

These laws express the conservation of mass, momentum and energy, 
respectively]]] where the vector (ijJo) is 

h ■= Wpj - W JP , l := Wp, t - W >tP . 

We can see that the momentum and energy are indeed conserved quan- 
tities: Iq and lj average to zero, since f(v* dp)p = f(v * p)dp for an 
even v. 

One can derive one more identity (a structure equation) by multi- 
plying the evolution equation by and taking the real part: 

A+ (- Ap + Wp) =0 . (11) 

This equation is the result of an infinitesimal transformation on the 
target when we regard as a map into the complex plane. Using 

1 As it is well known, any Euler-Lagrange equation derived from a local La- 
grangian has a conserved energy-momentum tensor; see, for instance, section 37.2 
in [4|. In our case, Tjk = (Jjk — &jk\ Tjq = —pj, Toj = <Joj and Too = — e. The 
vectors lo, lj are corrections due to the fact that our Lagrangian is nonlocal. 
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the structure equation, we can recast the conservation of momentum, 
equation fllObl) . into the form 

d tPj - V k {a/ + 5/( -Ap + Wp)} + lj = Q. 

Let us return to conservation laws (fIU|) . The conformal identity can 
be derived by contracting the mass equation ( 110 aft with |x| 2 /2; the 
momentum equation fllObl) with tx^\ and the energy equation (11 Pel) 
with t 2 ; and adding the resulting identities. The final result can be 
written in the abstract form 

d t e c - VjT j + r = , (12) 

where the relevant quantities are 

e c := (\x\ 2 /2)p + tx j p j + t 2 e = t 2 (J V (e i|a;|2/4 V) | 2 + Wp) , 

T j := {\x\ 2 /2)p> + tx k a\ + tx j (-Ap + Wp)+ t 2 a 3 , 
r := t\ + tx j lj - ntAp + t{n - 2)W p , 
By integrating (Tl2|) in space, we obtain the ODE 

E c + R c = , 

where 

E c := J dx{e c } , (13) 

E. := j ix{{ n-2)W P ^l i} ; (14) 

note that E c is the pseudoconformal energy. Next, we recast R c into a 
convenient form. By inspection of fTl4|) . it remains to compute 

l dx {Wx j pj - x j Wjp] = 2 / dx\dx 2 {f 1-2(^1 • Vi_ 2 (piP2)} 

= y tfoida^ { ( - 27Wi_ 2 - 2(x! - a; 2 ) • Vwi_ 2 ) (P1P2) } • 

In the above calculation, we used the fact that vi- 2 is symmetric 
with respect to the 1—^2 and 2 — > 1 transposition, while Vi_2 is 
antisymmetric. Moreover, we have 

Vi_2 ■ {xi - x 2 ) = 2n ; (x 1 - x 2 ) ■ Vi_ 2 ^i-2 = 2(x x - x 2 ) ■ Vf i_ 2 • 
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Finally, regarding ( fl4l) . notice that 

J dx {(n - 2)Wp} = J dx x dx 2 {2(n - 2)v 1 ^ 2 p 1 p 2 } • 
Substituting back into ffl4|) . we wind up with the integral 



R c = tj dx x dx 2 { [(-4)ui_ 2 - 2r 1 _ 2 v[_ 2 ] (P1P2)} • (15) 

This integral is used as an alternate expression for R c . 
Thus, we have proved the following lemma. 

Lemma 2.1. Let <fi be a solution of the Hartree equation (151) , and let 
E c , R c be defined by (113p . (|T5l) . Then, the following equation holds: 

E c + R c = . (16) 

Remark 2.2. In order to obtain a decreasing pseudoconformal energy, 
we need R c > 0, which is unfortunately not true for the Coulomb 
potential. Instead, we proceed to show that R c is integrable in time. 

We first state another consequence of our previous calculations. 

Lemma 2.3. Define 

E cc := f = jdx{t{\V (e<N 2 / 4 V) | 2 + W p) } . 
Then, E cc satisfies 

E cc + R cc = , 

where R cc is defined by 

Rcc := I V (e^l 2 / 4 V) | 2 -2 J (v + r^ 2 v' x _ 2 ) Pl p 2 . 

Remark 2.4. Notice that for v{x) = x(\ x \)m, Rcc is positive if xij) is 
decreasing for r > 0; thus, E cc is decreasing. 

In conclusion, using the Sobolev embedding and interpolation we 
have the following corollary. 

Corollary 2.5. Let <fi be a solution of the Hartree equation ([8]). Then, 
the following estimates hold for all t > 1 : 

||0(V)IU«>(R3) < ^^cc(l) , 

\m,-)\\ L , m <^E cc (i) . (17) 
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Using Lemma (12. II) . the result of Corollary 12.51 can be improved: 
Theorem 2.6. Let <p be a global smooth solution of the Hartree equation 

+ A0 - (t; * |0| 2 )0 = (18) 

with initial condition <fio such that E c (l) is finite. Then, 

\m,.)\\ L6 <ct-^ , 
\m,-)\w<ct-^ ie , 

and, thus, 

poo 

\m,-)\\ 2 L6m dt<c , 



U(t r )\n, m dt < c . 

Proof. Using the fact that — Av — 2rv' G L 1 together with (fT7|) . we see 
that 

Rc(t)<CtU(t,-)t L ,<Ct- 1 ' 2 . 
By integrating ([15]) . we conclude that 

E c (t) < Ct 1/2 

for t > 1. Using the Sobolev inequality and the definition of E c (see 
(|T3|) ) we conclude that 

||0(t,-)IUa<C||v(e^ 2 / 4 V) \\ L ,<Cr^ . 

Interpolation with energy conservation gives 

\m,-)\\L*<cr 9 / 16 . 

□ 



3. A priori Estimates 

In this section, by using Theorem 12.61 we derive a priori estimates for 
the solution u of (jlj). We recall definitions ([5]) of Theorem II .![ suitably 
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abbreviated, and add a few new ones: 

m 12 := = (~v(x - y)^(x)^(y)) , 

012 := -Ai5i2 + w 12 , 

Wn := -Wl-20102 > 

sh(fc) := it , 

ch(/c) := 1 + p := 5i2 + p , 
r := , 
q := uu . 

These are all operators kernels, and their products are interpreted 
as compositions. Notice that w and m have the symmetries u^i = w^, 
i.e., w* = to; and m,2\ = mi2, i.e., m T = m. The evolution equation for 
u = sh(fc), given by (jl]), is abbreviated to 

S(u) — (1 + p)m = + urn) ru , (19) 

where 

S(u) := m t + (?M + M(yf T , 

W(p) := ip t + [g,p] , 

and uu is symmetric, M21 = W12, i-e., u T = u, while P12 is self-adjoint, 

P21 = P12, i-e., P* = P- 

Notice that q is related to p by 

q = 2p + p 2 . (20) 
Trigonometric identities such as (120]) follow from e K e~ K = I for 

The key observation in this section is the following lemma. 
Lemma 3.1. The following identities hold: 

(W(p) + urn) r + r (W(p) — mu) = , (21) 

F := W(q) = mu(l +p) - (1 + p)um . (22) 

These equations are equivalent for any positive semi-definite kernel p, 
q = p 2 + 2p, and r = (1 + p)~ l . 

Proof. To prove (j2"T|) . multiply ffl9l on the right by u, take the adjoint 
of ffT9]) , namely, 

S'(w) — m(l + p) = wr (— W(p) + mu) 
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multiply it on the left by u, and then subtract. The resulting equation 
reads 

W(q) = W(p)rq + qrW(p) + umrq — qrrriu 
+ (1 + p)mu — um(l + p) . 

Now we can combine two terms as follows, using the fact that r = 
(1 + pY 1 and q = (1 + p) 2 - 1: 

umrq — «m(l + p) = um ((1 + p)~ [(1 + p) 2 — l] — (1 + p)) = —umr . 
Using the equation q = 2p + p 2 , we obtain 

W(2p + p 2 ) + umr — rrriu = W(p)r(2p + p 2 ) + (2p + p 2 )rW(p) ; 
hence, 

2W(p) + pW(p) + W(p)p + umr — rmu 

= W(p)(l + pY 1 [1 + (1 + p)]p + p [1 + (1 + p)] (1 + p)- 1 ^^) . 

Thus, we have 

2W(p) + umr — rmu = W(p)rp + prW(p) . 

Now observe that 1 — rp = r to conclude the proof of ( |2TT) . To prove 
( 1221) . multiply ( 12TT) on the right and left by (1 + p) and recall q = 
2p + p 2 . □ 

We are ready to state and prove our main a priori estimate: 

Theorem 3.2. Let u = sh(k) be a solution of fl4]) on some interval 
[0,T]. Then, the following estimate holds: 

\HT)\\ L2 < ( I \\m\\ L 2dt + \\u{Q)\\A exp I / \\m\\ L *dt I . (23) 




Proof. Taking the trace in (122]) we obtain 

_ \\u\\\i = tr [(l/i)(mu(l +p) - (1 +p)ufn)] . (24) 
Thus, we have 

jllMli 2 - 2 (ll m IU 2 ll M IU 2 + ll m IU 2 ll M IU 2 lblU 2 ) 

< 2 (||m||i2||u|| L 2 + ||m|| L 2 ||«||^a) • 

The inequality < follows by talking the trace of (120]) together 
with the observation that tr(p) > 0. Now we can employ a Gronwall 
type inequality to deduce ( 123]) . 

□ 
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Summarizing the results of the previous two sections, we draw our 
main conclusion. 

Corollary 3.3. Let <fi be a solution of the Hartree equation satisfying 
the assumptions of Theorem \2.6\ and let u = sh(k) be a solution of 
equation fl4]) on [0, T], as in Theorem \3.2[ Assume the potential v is in 
L 2 (R 3 ). Then the following estimate holds: 

\\u(T)\\ L i<C(l + \\u(0)\\ L 2) . (25) 

Proof. By Theorem 13.21 it suffices to control \\m\\ L i( dt ) L 2( dxdy y Notice 
that |H| 2 2(R6) = / (v 2 * \<p\ 2 ) \<p\ 2 dx < Cll^llii^ll^H^^g). Using the 
estimates of Theorem 12.61 we conclude that \\m\\ 11^12^^ < C. 

□ 

4. The local existence Theorem for equation (J4]) 

In paper I, we showed that (jl]) has local solutions provided u(0) = 
and v(x) = for \ Co°- 111 this section we relax the assumptions 
to u(0) G L 2 (IR 6 ) and v(x) = and prove local existence in an 

interval where W^W^do^DL^idxdy) is small. Notice that by Theorem I2.6[ 
[0, 00) can be divided into finitely many such intervals. This implies 
global existence for equation (j4j). 

In this setting, we can no longer assume that ||w||i,°°£,2 is small, and 
terms such as W{p)ru are no longer small compared to S{k) (see (TL9l) 
for the notation). Our equation seems quasilinear, but can be trans- 
formed into a semilinear one. In order to prove local existence, we 
must solve for u = sh(k) rather than k, and express p = y/1 + uu in 
the operator sense. Thus we have to prove the following proposition: 

Proposition 4.1. The map 

k 1 — y sh(k) = u 

is one to one, onto, continuous, with a continuous inverse, from sym- 
metric Hilbert- Schmidt kernels k onto symmetric Hilbert- Schmidt ker- 
nels u. 

Proof. The appropriate context for this proof is set by noticing that 
the equation u = sh(k) is equivalent to 

exp (° *\ = (^+^ «) 

^ \k OJ \ U y/l + Uu) 

By the spectral theorem, the exponential map is a continuous bijec- 
tion from self-adjoint Hilbert-Schmidt "matrices" to positive definite 
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"matrices" P for which \\I — P\\l 2 is finite. Our target matrix is 

p _ f Vl + uu u 

\ u Vl + uu 

Besides being positive definite, this matrix is symplectic; thus, it sat- 
isfies P T JP = J where 

-Co 1 ). 

and also satisfies LPL = P^ 1 where 

L-( l ° 
L ~ \U -1 

Thus, we have P = e p where p is self-adjoint. Since eP Je p = J, we 
conclude that p is symplectic, or p T J+Jp = 0. (Proof: e JpJ J = Je~ p J 
is always true; thus, by easy algebra e p = e JpJ . Since both p T and 
JpJ are self-adjoint, the exponential is one-to-one, and we conclude 
that p T J + Jp = 0.) Similarly, from Le p L = e~ p we infer LpL = —p. 
The first two conditions force p to be of the form 

fa b 

where a = a* , b = b T ', c = b* . The third condition entails a = 0. Thus, 
p can be re-written as 



P 



k 
k 



□ 

The main new ingredient of this section is the following theorem. 

Theorem 4.2. The following equations are equivalent for a symmetric, 
Hilbert- Schmidt u: 

S{u) = (l+p)m+ (W(p) + um)ru ; (26) 

S(u) = (1 + p)m + -[W(p), r\u + - (rmu + ufnr) u ; (27) 

S{u) = (1 + p)m + - [W, r\u + - (rrriu + umr) u ; (28) 
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where we set 

F := mu(l + p) — (1 + p)um , 



W:=— [ (q-z) l F(q-z) 1 VTTzdz, 
2iri J v ' v ' 



q := uu 



1 + p := + «« , 
r := (1 + ml) -1 . 

i/ere, r is a contour enclosing the spectrum of the non-negative Hilbert- 
Schmidt operator uu. Equation (|26|) zs the same as (jlj), suitably re- 
written. Note that F corresponds to W(q) and W corresponds to W(p). 

Proof. Assume u satisfies ( |26|) . Recalling the estimate (I2T|) we conclude 
u satisfies 

S(u) = ^ (W(p)r + rW(p)) u + - (W(p)r - rW(p)) u 
+ (1 + p)m + umru 

= - (rmu — ufnr) u + r\u + (1 + p)m + umru 

= -[W (p) , r]u H — (rmu + ufnr) u + (1 + p)m . 

2 2 

Thus, m satisfies (|27|) . Notice that both (B7I) and are of the form 

5(tt)=XM+(l+p)m, (29) 

where X is self-adjoint. To see that X is self-adjoint, notice that both 
W(p) and are skew-Hermitian. Then, the procedure can be reversed 
to show that if u satisfies (|2"9"|) then the identity ( T2"2"|) . and thus f l2T|) . 
holds. 

Indeed, composing the complex conjugate of fl2T|) on the left with u, 
we obtain 



uS(u) = uXu + u(l + p)m 
The adjoint of this operator is 



S(u)u = uXu + m(l + p)u 
Subtracting the first equation from the second one gives 



W(uu) = m(l + p)u — u(l + p)m 



which is the same as ( |22|) . using (1 + p)u = u(l + p) and u(l + p) - 
(1 +p)u. Thus, ( 126|) and (|27|) are equivalent, and all three equations 
(I26]h (127|) and (12^1) - imply the equivalent formulas (|2il, (l22|). 
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Next, assume ( 1271) holds. Then, we have [22] 

q = / (z — q)^ 1 dz and 

27TZ J 

r 

= f(q - z)- l VTTzdz , 

2m J 

r 

and 

- z)" 1 ) = -(g - zY^W^q - , 

^(V^T^) = -^ l{q-zY 1 W(q)(q-zY 1 VT+~zdz . (30) 
27T2 y 

r 

So, (J2HD follows, since W(p) = W(^TTq) and = F. 

Conversely, assume (128]) holds. Then, W(g) = F as before; and 
W(p) is given by Q3Qj), thus fl27J) holds. 

□ 



Theorem 4.3. Using the same notation as in Theorem 4-2, let uq G 



L (M b ) 6e symmetric, given. There exists e suc/i i/iai «/||^||Li([o,T])L 2 ((ixd3/) 
eo t/ien there exists u G L°°([0, T])L 2 (dxdy) solving ( |28l) with prescribed 
initial condition u(0, x, y) = Uq(x, y) G L 2 (R 6 ). TTie solution u satisfies 
the following additional properties: 

(1) 

i-Q^ — A x — A y I w||il([o,T])L 2 (efa(ij/) < C ! (31) 



< 



(2) 



.9 \ 

i-QL - A x + A y j p||ii([o,T])L 2 ((fa(iy) < C • (32) 



In i/iis context, p is defined as \/l + uu — 1. 

Proof. The equation (1281) is of the form 

S(u) =m + N(u) , (33) 

where N(u) involves no derivatives of u. Recall the fixed time esti- 
mate ||/c/||l 2 < II k || o P \\l || l 2 ) where op stands for the operator norm, 
and L? stands for the Hilbert-Schmidt norm. Since r and (q — z)~ x 
have uniformly bounded operator norms and ||p||l 2 < IMU 2 > an d also 
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\z\ < C||u||^2 on T, we have 

\\N( U )\\ L i L 2 < C(l + IklH^)!^!!^ , 

\\N(u) - N(v)\\ L i L 2 < Cmax{l, ||m||1oo L2 , IMIl°o L 2}||m|| L i L2 , 

X \\u — v\\loc, L 2 

where L x l? stands for ^([0, T])L 2 (R 6 ) and L°°L 2 = L°°([0, T})L 2 (R 6 ). 
Recalling the energy estimate 

IMIl°°l 2 < \\u(0, -)||l 2 + \\Su\\ L i L 2 , 

we see that, for any given C there exists an eo such that (|33|) has a 
fixed point solution in the set ||w||l°ol2 < C provided ||m||ii/_,2 < eo- 

To prove (I3~TI) . we already know that ||S'm||lil2 < C, so we must only 
account for the lower order terms in g, namely V\2$\$iU (composition 
of kernels) and (v*\<j)\ 2 )u (multiplication). These are both easy because 
we know ||«||l°°l2 < C and Theorem 12.61 implies ||ui20i02||i l x 2 < C* as 
well as ||f * 1 0| 2 ||z,ii°o < C , since v G L 2 . 

A similar proof applies in order to show that (1321) follows from esti- 
mate (22} . 

□ 

5. Estimates for error terms 

In this section we obtain estimates for the error terms J Q T || e^X^e - " 8 ^! || j dt 

(quartic term) and J Q T ||e B [A, V^]e — ^Ojl jr dt (cubic term). These terms 
were encountered in paper I. 

We start by recalling the following result (Proposition 2, section 7 of 
paper I): 

Proposition 5.1. The state e B Ve~ B Q has entries on the zeroth, sec- 
ond and fourth slot of a Fock space vector of the form given in paper I. 
In addition, if 

II ( — A x — A y I u\\ L i[o : T]L 2 (dxdy) < Cl, 

|| - &x + Ay^j p\\Li[0,T}L 2 (dxdy) < C 2 

and v(x) = x( x )^\) orv(x) = m ; then 

[ \\e B Ve- B tt\\ 2 T dt < C , 
Jo 

where C only depends on C\ and C 2 ■ 
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Based on this result, estimates ( )3T|) and ( 1321) and Cauchy-Schwarz 
in time we conclude: 

Proposition 5.2. The following estimate holds: 

f \\e B Ve- B n\\ T dt < CT 1 / 2 . (34) 
Jo 

Now we turn attention to J Q T ||e B [yl, V}e~ B £l\\jrdt, seizing the oppor- 
tunity of improving on results in section 8 of paper I. There, we had to 
estimate a certain trace; see equations (61) and (62) of paper I. This 
task can be avoided by commuting a X2 and a* 2 in equation (60) of paper 
I. Thus, terms involving sh(k)(x, x), as in (62) of paper I, can in fact be 
avoided. To illustrate this point, we include the calculations here in the 
Appendix, which in effect replaces section 8 of paper I, incorporating 
the above remark. Our result is now simpler and stronger. 

Proposition 5.3. The state e B [A,V]e~ B Q has entries in the first and 
third slot of a Fock space vector of the form ipj - ipjjj and ipj' - ipni' 
given in the Appendix. In addition, if 

f.d \ 

II — A z — Ay j u\\ L i^ T } L 2( dxdy } < Ci, 

f.d \ 

|| — A x + A y J p||ii[o,T]i 2 (da;d2/) < C 2 , 

II ^— + (pW^iO^L^idxdy) < C3 , (35) 

and v(x) = or v{x) = w ; then we have 

[ \\e B [A,V}e- B Q\\ 2 c dt < C , 
Jo 

where C only depends on G\, C 2 and C 2 . Thus, the following estimate 
holds: 

f \\e B [A, V]e~ B n\\ T dt < CT 1 / 2 . (36) 
Jo 

Remark 5.4. Notice that (135 p is satisfied by Theorem ( 12. 6p . 



Estimates ( 134|) and (|36|) form the basis of Theorem II. 3[ which is the 
main result of this paper. 
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6. The operator e B 
In paper I, we used the definition 

B(t) := ^ J (k(t,x,y)a x a y -k(t,x,y)a* x a* y ) dxdy (37) 

with ||/c||L 2 (da;(ij/) small; e B was defined as a convergent Taylor series on 
the dense subset of vectors in J 7 with finitely many nonzero compo- 
nents, and then it was extended to J 7 as a unitary operator. Consider 
the Lie algebra sp(R), or sp(C) of symplectic matrices with real (or 
complex), bounded operator coefficients. These satisfy JS + S T J = 
and have the form 



S 



a b 

t 



_ c —a 

where b = b T , c = c T . Further, consider the corresponding groups 
Sp(M), Sp(C) of bounded operators G which satisfy G T JG = J. In 
applications, G = e s G Sp is defined by a convergent Taylor series. By 

defimtion, G acts + * by actiug on the vector g) and, of 

course, preserves the symplectic form ^ J 4>ip. 

The following Lie algebra isomorphism from sp(C) to operators (not 
necessarily skew-Hermitian) was a crucial ingredient in paper I: 



d k \ T f d k 

I - d T}^- L \l _ d T 



):=\{a x o;)(f _ k dT )j(<) (38) 

/a x a* + a*a x If 
d{x,y) y —^— - — dxdy + - / k(x,y)a x a y dxdy 

~ \ J K x ,y) a l a * y dxd V ■ 

To ensure that the resulting operator is skew-Hermitian we now restrict 
this isomorphism to the Lie subalgebra sp c (R) := Csp(M)C T for 

This is a change of basis that will be explained below. 
Lemma 6.1. The map 

S ^ CSC T (39) 
is a Lie algebra isomorphism of sp(C) to sp(C). 
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Proof. The "matrix" C is unitary (C T = (C) _1 ) and also satisfies 
C T JC = CJC T = iJ; thus, i~%C and i^^C T belong to the symplectic 
group Sp(C) (and the choice of does not matter). Since (l39il does 
not change if we replace C by i~zC, we see that fl39|) is just conjugation 
by an element of Sp(C), and thus is a Lie algebra isomorphism. □ 

Lemma 6.2. If S G sp c (R) ; then I(S) is skew-Hermitian. 

Proof. This proof will also motivate the choice of C. 
Define the self-adjoint operators of "momentum" 

Px .- Dx =- ^L±< 



V2 



and "position" 



Q. := X x = ' 



✓2 

These satisfy the canonical relations 

[D x ,X,] = -5(x-t/) . 

We will rewrite fl38|) in terms of the self-adjoint operators .D and X. 
The change-of-basis formula is 

" A c(y') (40) 



a* ; I X, 

for 

i A -2 



see page 174, (4.13) of [11] for a closely related construction. Notice 
that JC = iCJ with C = (C T )~ 1 ; thus, 



2 

" 2 



At this point it is natural to introduce the Lie algebra isomorphism 

sp(C) ->■ sp(C), 

A = A a , a * -> A AJC := C T A a , a »a . (41) 
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Since C is unitary (C T = (C) -1 ), this is the inverse of ( )39|) . At this 
stage it is clear that if 

is real then the corresponding operator is skew- Hermit ian. Thus, sp c (R) 
consists of those A a>a * such that the corresponding A D X £ sp(M). □ 



Remark 6.3. In particular, for our K, 

K 

the corresponding decomposition in D x and X x is (see (HOI) ) 

%k -m 

thus, K G sp c (M). 

It is easy to check that, if S G sp c (R), then 

e 5 ( , ] is of the form ' ' 



k 
k 



Jf A * = C^ITC = (J* : (12) 



Thus, it is legitimate to parametrize the vector by and denote 

e 5 (0) := e s 

We now recall the results of section 4 in paper I: 

Theorem 6.4. Let (ft e L 2 and R, S G sp(C) unt/i L 2 for Hilbert- 
Schmidt) coefficients. Then 

[1(S),A^] = A(S^) , (43) 

and therefore 

[I(S),A(<f>)) = A(S(<t>)) , (44) 

[X(5),X(i?)]=X[5,i?] . (45) 
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In addition, if S G spc(M) and \\S\\l2 is 



e 1 



e~™ = A(e s ( f \ , (46) 



e^/(i?)e- I ( s )=X(e 5 i?e- s ) 



(47) 
(48) 

(|^)e-^-x((|^)e-). (49) 

Proof. The formulas (I43l)-( I45|) are elementary calculations. Formulas 
(I46l)-(l49l) follow by analyticity (power series) since e^ 5 -* is given by 
a convergent Taylor series on the dense subset of Fock space vectors 
with finitely many non-zero components. Replace S by tS (t G C, 
small) and check that all derivatives of the left-hand side agree with 
all derivatives of the right-hand side at t = 0. □ 

For 1 1 | £2 ( rf:r dy ) large, S G sp c (R), the series defining e x ^ may not 
converge on a dense subset. So, we define 



e 



AS) _ (JC(S)/n 



where n is so large that e x ^ s ^ n is defined by a convergent series on 
vectors with finitely many components, and is then extended as a uni- 
tary operator to J 7 . This definition is clearly independent of n and still 
satisfies the crucial properties (I4l)]) - (j4"9"]) . 

For the rest of this section, we discuss connections with well-known 
results and explain the change-of-basis formula. 

6.5. Connection to the Heisenberg group and metaplectic rep- 
resentation. Recall that the classical Heisenberg group H n is C n x 
M. with multiplication law (z,t)(w,s) = (z + w,t + s — Qzw); see 
(1.20) in [11] • In our setting, H is L 2 (M 3 ) x R with multiplication 

law (0, t)(ip, s) = (<j) + t/j,t + s - ^(jnp). The map ((f), t) ->■ e ~ AW e i4 is 
a unitary representation of H. Indeed, we have 

e -A{4>) e it e -AW) e is _ e -A{<t>+i>)+l[A(<t>),A{il>))jt+is 

_ e -A{<t>+Tl>)+\ J (Hi-Hj) e it+is 

_ e -A{4>+V) e i{t+s-$s5 <H>) _ 

Shale [21] extended the standard construction of the metaplectic rep- 
resentation (see chapter 4 in (TTJ) to the infinite dimensional "restricted 
symplectic group" rSp(R) = {T G Sp(R), (T*T) 1 / 2 -/is Hilbert-Schmidt}. 
We do not use his results directly; and the following comments are just 
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meant for completeness. His results, and those of [TT], are written with 
respect to the basis D X ,X X . By assuming G € rSp{R), Shale showed 
there exists a unitary transformation of J 7 , Y(G), such that 

e -A(G4>) = Y(G)e- AW Y(G)- 1 ; (50) 

also, any two such unitary transformations Yi(G), Y^iG) are related 
by Yi{G) = e Le Y2{G). The mapping G h-> Y(G) is a projective unitary 
representation, meaning that Y(Gi)Y(G 2 ) = e i6( - Gl ' G ^Y(GiG 2 ). In 
particular, we identify our unitary operator e x ^ (after we reconcile 
the bases) as e x ^ = e t9 Y(e s ) for some 9 = 9(S) E R; we skip further 
details. 



Appendix A. Computation of cubic error term 

With recourse to equation (56) of paper I, and because of the com- 
ments following Proposition (15.21) . we now carefully compute the error 
term 

e B [A, V]e- B = J v(x - y) (^>(y)e B a x e~ B e B a x e~ B e B a y e~ B (51) 

+ (fi(y)e B a* x e~ B e B a*e~ B e B a x e~ B J dx dy , 

(52) 

which acts on the vacuum state, Q. All terms ending in a can be 
ignored. After commuting all a terms to the right, we are left with a 
pure cubic and a pure linear term in a*, which we proceed to compute. 
Recall the following formula proved in paper I: 

(53) 
Thus, we have 

and, similarly, 

e B a* x e~ B = I \sh(k)(y, x)a y + ch(k)(y, x)a*) dy . 



e [a y ,a y ) 



-B 



[a y ,a y )e' 
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We are ready to extract the pure a* term from ( l5TT) . Before any 
simplifications, (JoTj) reads 



v(x-y)(f)(y) [sh(k)(z 1 ,x)a Zl + ch.(k)(z u x)a* zi 



ch(k)(z 2 + s\i{k)(z 2 , 



ch(k)(z 3 ,y)a Z3 + sh(k)(z 3 ,y)a* Zz ) dz 1 dz 2 dz 3 dxdy . 
Thus, ( 151]) contributes the cubic term 



/ = y v(x-y)(j)(y) \ ch(k)(z 1 ,x)sh(k)(z 2 ,x)sh(k)(z 3 ,y)Ja* zl a* Z2 a* Z3 

dz\dz 2 dz 3 dxdy . 

Thus /(^) has entries in the third slot of Fock space equal to (after 
normalization and symmetrization) 



ijji (zi,z 2 ,z 3 ) = J v(x-y)(f)(y) \ ch(k)(z 1 ,x)sh(k)(z 2 ,x)sh(k)(z 3 ,y)jdxdy . 

For the linear terms, keep only the a* term from the last row, and 
exactly one a and one a*s from the first and second rows, and commute 
the a's to the right. Hence, we are left with two terms: 



II = I v(x — y)(f)(y)sh(k)(zi, x)sh(k)(z 2 , x)sh(k)(z 3} y)a Zl a* Z2 a* z dz\dz 2 dz 3 dxdy 



v(x - y)(p(y)sh(k)(z 1 ,x)sh(k)(z 2 ,x)sh(k)(z 3 ,y) 
(S(zi — z 3 )a* Z2 + 5{zi — z 2 )a* Z3 j dz\dz 2 dz 3 dxdy (modulo linear terms in a) 



v(x — y)(f)(y)sh(k)(zi, x)sh(k)(z 2 , x)sh(k)(zi 1 y)a* Z2 dzidz 2 dxdy 



+ / v(x — y)(f)(y)sh.(k)(zi,x)sh(k)(zi,x)sh(k)(z 3 ,y)a* z . i dzidz 3 dxdy 

and 



/// = / v(x — y)(j)(y)ch(k)(zi,x)ch(k)(z 2} x)sh(k)(z 3l y)a* Zi a Z2 a* Z3 dzidz 2 dz 3 dxdy 



v(x — y)(J)(y)ch(k)(zi } x)ch(k)(z 2j x)sh(k)(z 2 , y)a* zi dz\dz 2 dxdy . 
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These terms contribute to the first slot of Fock space entries of the 
form 



ijjjj(z) = / v(x — y)(f)(y)sh(k)(z 1 , x)sh(k)(z, x)sh(k)(zi,y)dzidxdy 



+ / v(x — y)(j)(y)sh(k)(zi,x)sh.(k)(zi,x)sh(k)(z, y)dz\dxdy 

and 



ipin(z) = = J v(x — y)4>(y)ch.(k)(z,x)ch.(k)(z 2 ,x)sh(k)(z 2 ,y)dz 2 dxdy . 
Next, we concentrate on the contributions of (1521) . 



52]) = / v(x -y)(j){y) [sh(k)(z 1 ,x)a zi + ch(k)(z 1 ,x)a* Zi 



sh(k)(z 2: y)a Z2 + ch(k)(z 2 , y)a* z 



■-■2 



ch(k)(z 3 ,x)a Z3 + sh(k)(z 3 ,x)a* Z3 J dz\dz 2 dz 3 dxdy 
which contributes a cubic in a*: 



I' = J v(x — y)(f)(y)ch(k)(zi,x)ch(k)(z 2 ,y)sh(k)(z3,x)al i al 2 a* Z3 dzidz 2 dz3dxdy . 
Here, we have 



ipriz^z^zs) = J v(x - y)(f)(y)ch(k)(z 1 ,x)ch(k)(z 2 ,y)sh(k)(z 3 ,x)dxdy 
The linear terms in a* (modulo a) are 



//' = / v(x — y)<p(y)sh(k)(zi, x)a Zl ch(k)(z 2 ,y)a* z sh(k)(z 3 , x)a* dzidz 2 dz 3 dxdy 



v(x — y)4>(y)sh(k)(zi, x)ch(k)(z 2 , y)sh(k)(zi, x)a* Z2 dz\dz 2 dxdy 



v(x — y)(f)(y)sh(k)(z 1 , x)ch(k)(zi, y)sh(k)(z 3 , x)a* z dz\dz 3 dxdy 



and 



III' = / v(x — y)(j)(y)ch(k)(zi,x)sh(k)(z 2 ,y)sh(k)(z 3 ,x)a* Zi a Z2 a* Z3 dzidz 2 dz 3 dxdy 



v(x — y)(J)(y)ch.(k)(zi, x)sh(k)(z 2 , y)sh(k)(z 2} x)a* Zl dzidz 2 dxdy . 
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The Fock space entries read 

ipir(z) — J v(x — y)<f)(y)sh(k)(zi,x)ch(k)(z,y)sh(k)(zi,x) dzidxdy 



+ J v (x — y)(f)(y)sh(k)(zi, x)ch(k)(zi, y)sh(k)(z, x) dz\dxdy 

and 



^iii'(z) = / v(x — y)(f)(y)ch.(k)(z,x)sh(k)(z2,y)sh.(k)(z2,x) dz 2 dxdy . 



All the resulting ip can be estimated in L 2 (dtdz) by the method of 
section 7 of paper I, without using X s ^ spaces. We remind the reader 
how to estimate these terms. Take, for instance, ipnr(t, z). Write 
ch(k)(t, z, x) = 5(z — x) + p(t,z, x) to express tynv — i>s + Vv We 
estimate the first of these terms: 



\M Z )\ = \ / v(z -y)(f)(y)sh(k)(z2,y)sh(k)(z2,z)dz 2 dy\ 



< \\v(z - y)(f)(y)sh(k)(z 2 , z)\\ L 2 {dZ2dy) \\sh(k)(z2,y)\\ L ^dz2dy) ■ 
The second term is uniformly bounded in time; thus, 



3 

oc 



ips{t,z)\ dtdz 

J \v(z - y)(j)(y)sh(k)(z2, z)\ 2 dzdz 2 dydt < C 

by a local smoothing type result (see Lemma 2, section 7 of paper I). 
All other terms can he estimated by the same method. 
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